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The MBE grown Sii_a:Gea; islands on reconstructed high index surfaces, such as, Si(5 5 12), 
Si(5 5 7) and Si(5 5 3) show a universahty in the shape evaluation and the growth exponent param- 
eters, irrespective of the substrate orientations and size of the island structures. This phenomena 
has been explained by incorporating a deviation parameter (e) to the surface barrier term {Ed) in 
the kinematic Monte Carlo (kMC) simulations as one of the plausible mechanisms. 



PACS numbers: 81.15.Hi, 68.35.-p, 68.37.-d, 66.30.Pa 

The implementation of Ge and Sii-a^Gca: nanostruc- 
tures into Si-based devices is of great potential for fu- 
ture high-speed devices, due to advantages like enhanced 
carrier mobility and smaller bandgap Growth 
of Ge islands on clean Si substrates (particularly on 
low-index-oriented) has been extensively studied as a 
model system to understand hetero-epitaxy and Stranski- 
Krastanov (SK) growth mechanism But, the study 

of Ge growth on high index surfaces, such as Si (5 5 12), 
Si(5 5 7) and Si(5 5 3) is an area where very limited re- 
search work has been carried out 0] . Previous study 
by Kim et al Q reported on relatively thinner Ge growth 
on Si (5 5 12) and did not observe any shape transfor- 
mation of the SiGe structures. For the Ge on Si sys- 
tem strain relief and diffusion play a role in determining 
the morphology and composition of Ge or SiGe struc- 
tures 0, 0, 0- the present work, we have observed 
the variation in the size of the Sii_a;Gea; structures de- 
pending on the substrate orientation while having the 
similar aspect ratios and growth exponent values for all 
three substrate orientations. The variation in the size 
of the islands has been explained by considering the ex- 
perimentally observed composition of the structures and 
strain associated in them. The universal nature of shape 
evaluation and exponent values are simulated with 2D 
kMC simulations. With the interesting reconstructions 
on high index surfaces and anisotropic diffusion domi- 
nating during the direct current heating conditions, the 
present work stimulates the self-assembly growth related 
work in the area of nanoscience and nanotechnology. 

Reconstructed high index silicon surfaces consisting of 
alternating terraces and atomic steps can be used as tem- 
plates to form aligned one dimensional (ID) nanostruc- 
tures (8l-[lo|. Atomic steps on the high index silicon sur- 
faces are responsible for many surface dynamic processes 
like surface migration and step diffusion. It is possible 
to use the high index substrates for the growth of self- 
organized nanostructures Q . Among the high index sili- 
con surfaces, Si(5 5 12) is oriented 30.5° away from (001) 
towards (111) with one-dimensional periodicity over a 
large unit cell width [l,[lH. Si(5 5 7) with vicinal angle 
of 9.45° from (111) towards (112 ) 12] and Si(5 5 3) , 



tilted at -12.27° from the (111) plane towards the (0 1) 
plane [13j, are important high index Si surfaces. In all 
the above vicinal surfaces, the step edges are along (110) 
direction. 

It is known that strained epitaxial layers tend initially 
to grow as dislocation-free islands and as they increase 
in size, may undergo a shape transition [tI. [i3]. Below 
a critical size, islands can have a compact symmetric 
shape. But at larger sizes, they adopt a long shape, 
which allows better elastic relaxation of the island's stress 
0j [3 EE] • Tersoff and Tromp proposed a model to ex- 
plain the growth kinetics involving a shape transition 
at a critical size by finding an expression for energy 
of dislocation- free strained islands [ij]. Also, Heyn re- 
ported the kinetic Monte Carlo (kMC) methods to study 
the influence of the anisotropy of surface diffusion and 
of the binding energies [l3]- Kinetic Monte Carlo meth- 
ods have been used to study the formation of nanoisland 
structures in a number of works [TtI . [18| , however, shape 
transition phenomena has hardly been studied. 

The experiments discussed in the following were per- 
formed in ultra-high vacuum(UHV), at a base pressure 
of ^ 2.5 X 10~^° mbar in a custom built molecular beam 
epitaxy (MBE) system [H. Samples of Si(5 5 12), Si(5 5 
7) and Si(5 5 3) were prepared from p-type boron doped 
wafers (of resistivity of 10 - 15 f2 cm). Substrates were de- 
gassed at 600°C for about 12 hours followed by repeated 
flashing (with direct current heating) for 30 sec. at a 
temperature of 1250°C to remove the native oxide layer 
to obtain a clean and well-reconstructed surface. The 
reconstruction has been confirmed with in-situ reflection 
high energy electron diffraction. The temperature was 
monitored with an infra-red pyrometer calibrated with a 
thermocouple attached to the sample holder. Ge was de- 
posited to various thicknesses of 2 to 10 monolayer (ML) 
at a typical deposition rate of 0.6 ML/min at substrate 
temperature 600°C [with direct current heating(DII)]. 
The samples were then post annealed at a temperature of 
600°C for 15 minutes by DH method. Also a set of sam- 
ples were prepared with radiative heating (RH) (where 
heating is achieved through a filament underneath). The 
post growth characterization of the samples was carried 



FIG. 3: SEM micrographs of (a)3 ML(b)5 ML(c)8 ML (d)10 
ML Ge/Si(5 5 3) at 600°C-DH 



FIG. 1: SEM micrographs of (a)2 ML(b)5 ML(c)8 ML (d)10 
ML Ge/Si(5 5 12) at 600°C-DH. Insets show magnified image 
of one structure. 




FIG. 2: SEM micrographs of (a)3 ML(b)5 ML(c)8 ML (d)10 
ML Ge/Si(5 5 7) at 600°C-DH 
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FIG. 4: Plot showing the Length Exponents of the strucrures 
formed on Si(5 5 12), Si(5 5 7)and Si(5 5 3) surfaces 



out ex-situ by field emission gun based scanning electron 
microscopy (FEGSEM) with 20 kV electrons. 

Figure [T] (a) and (b) depict the FEGSEM image for 
2 ML and 5ML thick Ge film deposited on Si(5 5 12) 
substrate. As shown in Fig (Ha), the island struc- 
tures are found to be spherical shape for a 2 ML Ge 
deposition, but for 5 ML Ge growth, rectangular is- 
land structures have been observed. For the 8ML and 
lOML deposition cases, we observe trapezoidal struc- 
tures [FiglTJc),(d)]. It was reported earlier that DH 
is a necessary condition to form trapezoidal structural 
structures on Si(5 5 12) The composition of these 

nanostructures has been characterized by using STEM 
- EDS 2^, 21 1, Rutherford backscattering spectrometry 



(RBS) 20] and the synchrotron-based high resolution x- 
ray diffraction (HRXRD). The HRXRD showed presence 
of graded Sii_j;Gej; system for the Si(5 5 12)and is also 
confirmed by RBS measurements. While RBS measure- 
ments show no prominent graded Sii-^Gcx structures for 
the case of Si(5 5 3)and Si(5 5 7) substrate orientations. 

In figlH Ge growth of overlayer 3ML to lOML on Si(5 5 
7) has been shown where spherical nanoislands and rect- 
angular nano rods are formed for 3ML and 5ML thick- 
nesses[Figl2Ua),(b)] and nano trapezoid in the case of 8 
ML and 10 ML Ge growth [FiglS^c), (d)]. We have main- 
tained the same growth condition for Si (5 5 3) substrate 
and seen that the growth kinetics and shape transforma- 
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TABLE I: Mean Length (L) 



Thickness 


Ge/Si(5 5 12) 


Ge/Si(5 5 7) 


Ge/Si(5 5 3) 


(ML) 


Mean L in (/im) Mean L in(nm) 


Mean L in(nm) 


3 


2.63 ± 0.14 


115.5 ± 7.2 


117.5 ± 8.2 


5 


3.51 ± 0.12 


200.6 ± 8.7 


189.8 ± 9.7 


8 


5.85 ± 0.22 


250.6 ± 9.6 


248.7 ± 11.6 


10 


6.25 ± 0.27 


297.8 ±11.4 


288.3 ± 12.4 


TABLE IL Aspect Ratios (length/ 


width) 


Thickness 


Ge/Si(5 5 12) 


Ge/Si(5 5 7) 


Ge/Si(5 5 3) 


(ML) 


Aspect Ratio 


Aspect Ratio 


Aspect Ratio 


3 


2.10 ± 0.1 


1.79 ± 0.2 


1.74 ± 0.6 


5 


2.20 ± 0.1 


2.28 ± 0.6 


2.34 ± 0.8 


8 


2.84 ± 0.5 


2.88 ± 0.3 


2.92 ± 0.9 


10 


3.13 ± 0.3 


3.17 ± 0.4 


3.15 ± 0.5 



tion follow in a similar way resulting from rectangular 
nano rod to nano trapezoid, as shown in figEl 

We have measured the length and aspect ratio 
(length/ width) of all the Ge-Si structures which has been 
shown in a tabular form in tableUand tablelTTl In the case 
of Ge/Si (5 5 12) system, the size of the aligned struc- 
tures are of micrometer size. But for the Si(5 5 7) and 
Si(5 5 3) systems, the size of the aligned structures are 
of nanometer size, though the shape evolution and also 
the aspect ratios are similar to that of Si(5 5 12) . From 
table U it is clearly seen that with increasing thickness, 
the size and aspect ratio of the Ge-Si structures increase 
accordingly. The aspect ratios of the Si-Ge elongated 
structures on Si( 5 5 12),Si(5 5 7) and Si(5 5 3) increase 
in a similar fashion as a function of increasing Ge growth 
coverages, which has been shown in table HIl 

The growth of island followed a universality and we 
found from our experiments that the mean length of is- 
lands {L) ^ t" where t is the time of deposition. The 
growth exponent from the experimental data is found to 
be about 0.75 ± 0.02 for the case of the direct heating 
and 0.5 ± 0.01 for the radiative heating case in all the 
three systems [table ITTT] .We observe that since a < 1, the 
growth of island is sublinear. In fig 21 the log-log plot 
for size(length) vs deposition time shows almost same 
slope(a) irrespective of the size of the Si-Ge structures 
on the three substrates. We now disscuss the shape trans- 
formation of structures by kMC. 

The kinetic Monte Carlo simulations were performed 
on a i X i square lattice with L = 100,200,300,400. 



TABLE IIL Length Exponent (a) 

Exponent Ge/Si(5 5 12) Ge/Si(5 5 7) Ge/Si(5 5 3) from kM C 
DH 0.77±0.04 0.76±0.03 0.73±0.02 0.60±0.02 
RH 0.51±0.03 0.53±0.02 0.52±0.04 0.45±0.02 



A coordinate system has been chosen such that the x- 
axis is directed along perpendicular to the step edges. 
We define horizontal bond and vertical bond as a pair 
of nearest neighbor atoms having the same y coordinates 
and same x coordinates respectively. Similarly a vertical 
bond is defined by a nearest neighbor pair of atoms with 
same x coordinates. The hopping rate of an adatom is 
given by 

where kB is the Boltzman constant, T is the tempera- 
ture and vq — iksT /h is the vibrational frequency in 
the direction of the hopping. A typical value for our 
case is vq = 3.6 x 10^^ sec~^. The quantities Ei, E2 
and i?3 denote the binding energies for the horizontal 
bond, vertical bond and the next nearest neighbor in- 
teraction respectively ; ni,n2 and are the number of 
horizontal bonds, vertical bonds and the number of next 
nearest neighbors respectively and Eb is a surface bar- 
rier term. For Eb and {i — 1,2,3) constant, we 
obtained spherically symmetric island structures in the 
kMC simulation, since the binding energies are same for 
all bonds and surface barrier term is independent of di- 
rection. To obtain asymmetrical structures , we need 
to break the hopping symmetry. In our model, we in- 
troduce anisotropy through binding energies of different 
types of bonds and the dependence of surface barrier on 
the direction of hopping. For the isotropic cases, how- 
ever, the hopping rate depends only on the nearest and 
next nearest neighbor interactions but does not depend 
on the specific arrangement of the neighboring atoms. In 
our experiment, the shape transformation has been ob- 
served when Ge is grown under DH condition [i^]. To 
simulate the experimentally observed shape variations, 
the surface barrier term is modified, so that it allows an 
asymmetric hopping of an adatom along the direction 
perpendicular to the step edge ( i.e. ± x direction). The 
values of En along the step edges are same in both ± 
y directions. When Ei ~ E2 and Ed does not depend 
on the direction of hopping, spherical islands are formed. 
Symmetric elongated structures are formed for E2 > Ei 
and Ed uniform in all the directions of hopping. Trape- 
zoidal islands can only be obtained when an asymmetric 
hopping term is present. The values of Ed are defined as 
follows. a) Ed = i?o(l + e) for hopping along the positive 
x-direction, b) Ed = Ea{\ — e) for hopping along the neg- 
ative x-direction, c) Ed = Eq for hopping along y-axis 
where e depends only on the magnitude of the current 
(but for the simulations, this is a parameter only). Also 
< e < 1 so that the probabilities are strictly less than 
unity. Note that when e = 0, it reduces to the case of 
anisotropic island formation for E2 > Ei as in the case 
of Ref . tl6i and isotropic islands formation for E2 = Ei as 
in Ref. [22| -Therefore, in our kMC simulation a nonzero 
value of e should correspond to the MBE growth done 
under DC heating. The mean shape of islands for this 
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FIG. 5: Snapshot of kMC island growth (a)spherical shaped 
islands at e = 0, £1 = E2 (b) e = 0.8 ,Ei < E2 rod like 
structures at coverage 0.26 (c) e — 0.8 , E\ < E2 showing 
shape transition by forming trapezoidal island structures at 
coverage 6. 



case is found to be trapezoids. 



We follow the kMC algorithm in [22,1231. We set the to- 
tal coverage 9 and the total number of Monte carlo steps 
A'mc- Particles are deposited at a constant flux 0/Nmc- 
We define a dimensionless scale parameter (p — Eq/KbT, 
where Eq is the surface energy barrier of the system when 
e = 0. This sets the time step for the simulation. An 
atom with empty adjacent site is called 'active'. An 
active atom is chosen at random and a single diffusion 
event is allowed to occur with a probability consistent 
with Eq.(l). If an atom is already there on the site to 
which it chooses to hop, then hopping fails. Time is incre- 
mented irrespective of whether the hopping is successful 
or not. The following parameters were used for our kMC 
: (j) = l.O^Ei = 1.0, £^2 = 6.0, E3 = 0.5, Eq = 1.0, e = 0.8, 
and coverage is 2.5%. It is important to note that the 
relative strength of the energy values is crucial in deter- 
mining the shape of the structures. The above set of 
energy parameters is one such example. 

In fig. [5] (a), spherical shaped islands for e — 0, Ei — E2 
are seen. For e = 0, Ei < E2 at coverage 0.29, rod like 
aligned structures formed [figlSJb)]. Figure|5] (c) is final 
snapshot of the trapezoidal structures formed at coverage 
for e = 0.8 and Ei < E2. From the kMC simulation, 
we have found a = 0.6 ± 0.02(for e = 0) and a — 0.45 ± 
0.03(for e = 0.8). 

In summary,we report the shape evolution of MBE 
grown Sii-^^Ge^: islands on reconstructed high index Si(5 
5 12), Si(5 5 7) and Si(5 5 3) surfaces. We show that a self 
assembled growth at optimum thickness leads to interest- 
ing shape transformations, namely, spherical islands to 
rectangular nanostructures and then to elongated trape- 
zoidal structures. We have experimentally observed an 
universality in the growth of the islands for all three high 
index surfaces by evaluating the aspect ratios and growth 
exponent. The growth exponent (for the longer side of 
the structures) experimentally found to be about 0.75 
± 0.02 for the case of the direct heating and 0.5 ± 0.01 
for the radiative heating case in all the three systems. 
Our kMC simulations show that such variations can be 
understood by introducing a deviation parameter e in 



surface barrier term E^. The experimentally observed 
shape variations and the growth exponent values are in 
good agreement with kMC simulations. This suggests 
the role of stochastic process involved in the shape tran- 
sitions of nanoscale structures. We thank S. D.Mohanty 
for stimulating discussions. PVS would like to thank the 
Department of Atomic Energy, Government of India for 
granting FEGSEM under 11th plan. 
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